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Abstract 

Using the Algebraic Bethe Ansatz we consider the correlation functions of the in- 
tegrable higher spin chains. We apply a method recently developed for the spin ^ 
Heisenberg chain, based on the solution of the quantum inverse problem. We con- 
struct a representation for the correlation functions on a finite chain for arbitrary 
spin. Then we show how the string solutions of the Bethe equations can be con- 
sidered in the framework of this approach in the thermodynamic limit. Finally, a 
multiple integral representation for the spin 1 zero temperature correlation functions 
is obtained in the thermodynamic limit. 
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Introduction 



A new method of computation of the correlation functions and of the form factors of the Heisen- 
berg spin ^ chains developed in ^ based on the Algebraic Bethe Ansatz [Q and the 

resolution of the quantum inverse problem [|l|, |5| has given a possibility to calculate a very 
large class of correlation functions for a finite chain and in the thermodynamic limit. The zero 
temperature correlation functions (which are defined as mean values of some products of local 
operators with respect to the ground state of the model) were obtained ^ as multiple integrals 
which coincide for the case without magnetic field with the results obtained using the vertex 
operators approach 0. However the new method gave a better understanding of the structure 
of these results. It was shown, in particular, that the expressions under integrals can be sepa- 
rated into two parts with different origin: algebraic part which depends on the choice of local 
operators and does not depend on the choice of the ground state and analytic (or determinant) 
part which, on the contrary, is fixed uniquely by the ground state. 

This very particular structure permits to hope that these results can be generalised for 
several more general situations in particular for the temperature dependent correlation functions 
for the Heisenberg spin ^ but also for other integrable models with the same iZ-matrix, and 
first of all for the higher spin chains (they can be considered as the first step to the future 
generalisations). These two apparently very different problems have however one common detail: 
the main difficulty is the analysis of excited states in the first case or of more complicated ground 
states for the second one. In both cases we should deal with bound states or, more precisely, with 
the string solutions of the Bethe equations. The understanding of the influence of these bound 
states (quantum breathers) is a very important step to the calculation of the temperature and 
time dependent correlation functions. 

For this reason, before considering a more complicated example of temperature dependent 
correlation functions, we consider the higher spin Heisenberg models. This case is more simple 
as the ground state of a higher spin XXX chain contains only strings of one kind and not a 
mixture of different types of strings as an arbitrary excited state. This problem is important 
also from other point of view as it can give some information about other integrable models 
including integrable quantum field theories. 

In this paper we consider the correlation functions of the XXX higher spin chains. This 
model was first considered by Kulish, Sklyanin and Reshetikhin in |p where the notion of fusion 
was introduced. It was solved by means of the Algebraic Bethe Ansatz by L. Takhtajan Q 
and independently by H. Babujian in the thermodynamics of this model was considered in 
jlOI I . The XXZ version of higher spin chains was introduced in [^] and solved in ||l^ , but these 
models are not considered here for several reasons. The correlation functions of the XXZ spin 1 
model in the anti-ferromagnetic regime were calculated in |13, 14, 15] using the vertex operator 
approach. We propose here a different way of calculation of the correlation functions based on 
the Algebraic Bethe Ansatz. 

As in the first step of computation of the correlation functions is the solution of the 
quantum inverse problem. Such a solution for a very large class of quantum integrable models 
including the higher spin Heisenberg chains was recently obtained Q in a form very similar to 
the spin ^ case. Taken together with the results for the scalar products of the Bethe states 
|16, 17, 18] it permits to obtain a representation for the finite chain correlation functions for 
arbitrary spin. On this stage one should take the thermodynamic limit and, hence, introduce 
the string solutions. We illustrate this procedure using the simplest example of the spin 1 chain. 
We would like to underline that a similar procedure is possible also for higher spin models but it 
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leads to more cumbersome calculations and we present only one form of the result for this case 
without detailed derivation. 

The main difficulty which arises from the presence of bound states is the fact that the 
algebraic part becomes singular. This problem can be solved by chosing the integration contours 
in the multiple integral representations taking into account the sign of the finite size corrections 



to the string picture |19, 20, 21 1. After this modification one can see that two parts appear again 
and the determinant part once again is defined uniquely by the ground state. Such a result is 
given in section 5 of this paper. 

The main result of this paper is the fact that the mean values with respect to the states 
containing bound states can be calculated in the framework of our approach. It means in 
particular that some new tools introduced here can be used also to calculate the temperature 
dependent correlation functions. 

This paper is organised as follows. In the first section we introduce the higher spin Heisenberg 
chains following the papers of Takhtajan |^ and Babujian [|lO|. The solution of the inverse 
problem for these models |^] is given in the second section. This solution is used to obtain 
representations for the correlation functions on a finite chain for arbitrary spin in section 3, we 
show, in particular, how to override the additional algebraic difficulties appearing in the higher 
spin case. The thermodynamic limit of the spin 1 chain is considered in two last sections. We 
show how to deal with 2-strings in the thermodynamic limit using the simplest example of one 
point functions in section 4. This first example permits two elaborate some simple rules which 
permit to deal with strings for general correlation functions (section 5). 

1 XXX Heisenberg chain with arbitrary spin 

In this section we introduce the XXX Heisenberg chains for arbitrary spin. We follow in general 
here the papers of Takhtajan and Babujian |p^ . 

Unlike the spin ^ case we start directly from the L-operator and construct later the Hamil- 
tonian from the transfer matrix. It is necessary to obtain an integrable generalisation of the 
usual Heisenberg chain (a direct generalisation of the spin ^ XXX Hamiltonian is not integrable 
for higher spins). However the L-operator can be obtained by the direct generalisation: 

^" A-^(.+i) V -^^^ A+z(.^+i); ^^-^^ 

One should note that for this L-operator the auxiliary space is two-dimensional but the quantum 
space has 2s + 1 dimensions. The matrices s^, are the spin operators in the representation 
of spin s. This L- operator has the same intertwining relation with the rational 4 x 4 i? matrix 
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(1.2) 



as in the spin ^ case: 

Ri2{\ - m) (L™)i(A) {LMii) = (L™)2(m) (L,„)i(A) i?i2(A - li). (1.3) 



2 



As usual, indices 1 and 2 in L operators denote two different auxiliary spaces. 

The next step is the construction of the monodromy matrix for a spin chain of M sites: 



r(A) = La/(A - ai)LM-i{X - ^M-i) . ■ ■ Li{X - ^i) = 



AiX) BiX) 
C{X) D{X) 



with arbitrary inhomogeneity parameters . For this monodromy matrix one can also write the 
"commutation relation" with i?-matrix: 

Ri2ix - ^) ri(A) T2{fi) = r2(/i) ri(A) i?i2(A - (i.4) 

and it contains the commutation relation for the operators ^(A), B{X), C(A) and D{X) acting 
in the quantum space. 

Up to this point all the construction was the same as in the spin | case but to construct local 
Hamiltonians from the monodromy matrix one should introduce some new concepts. The trace 
identities for the spin ^ case were based on the fact that Ln{0) is just a transposition operator 
in the tensor product of the auxiliary and quantum spaces. Now the quantum and auxiliary 
spaces have different numbers of dimensions. Hence it is necessary to construct a monodromy 
matrix with a 2s + 1 dimensional auxiliary space. Such a construction was proposed by Kulish, 
Sklyanin and Reshetikhin P] and it is called fusion. 

The fused L-operator L„'*(A) can be constructed as a following projection on the symmetric 
subspace in the tensor product of 2s auxiliary spaces (this symmetric subspace has exactly 2s + 1 
dimensions) of the following product of local L-operators: 

= Pi,...a„ia2.,n(A + 2is - i) . . . L„„„(A + i)La„n(A)P+,...,„ (1.5) 

where indices aj mark the corresponding auxiliary spaces and is the symmetrizer (pro- 

jector on the symmetric subspace). Thus we constructed the L-operator with 2s + 1-dimensional 
auxiliary space and satisfying the following property: 

where Pa^n is the transposition in the tensor product of the auxiliary space and local quantum 
space. This property is crucial for the construction of local Hamiltonians |1C] and for the solution 
of the inverse problem Q . The fused monodromy matrix is constructed as usual: 

T(^)(A) = 45 (A - - . • • l['\x - 6) 

The fused transfer matrix which is the trace of the monodromy matrix taken in the auxiliary 
space, 

TW(A)=tr„r„(A) 

commutes not only with the transfer matrices for any value of parameter A but also with the 
"ordinary" transfer matrix A{X) + D{X) (as it is a polynomial function of A{X) + D{X)): 

[r(^)(A),rW(/i)] = [tW(A), A(m) +D(^)] = 0. (1.6) 

It means, in particular, that the Hamiltonians constructed from the fused transfer matrix in the 
homogeneous case = using the trace identities: 



if(^) = const -^tW(A) 
dX 



(1-7) 

\=—is 
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also commute with the "ordinary" transfer matrix and can be diagonaUzed by the Algebraic Bethe 
Ansatz procedure. The Hamiltonians constructed by trace identities are local, translation 
invariant and can be written as polynomials of degree 2s of the local spin-spin interaction terms 
H: 



M 



H^'^ = Q2s{SmSm+l), (1.8) 

m=l 
2s / j .\ 2s 

j=l \k=l i-o ^' 

where s„ = (s^, Sn, s^) are spin operators and xi = + 1) — 2s{s + 1)]. The first example of 
such a Hamiltonian is the spin 1 case where 

-ff*-^^ = ^mSm+l - {SmSm+lf- (1-10) 

m=l 

As we already mentioned to diagonalize this Hamiltonian one can use the usual Algebraic 
Bethe Ansatz. We start from the ferromagnetic state | ) with all the spins up 

s+|0) = 0, Vn 



which is an eigenstate of the Hamiltonians (|l.7| ) and consider the action of the generalised 
creation operators B{X) on this state. It is easy to see that a state 

B{Xi)B{X2)...B{Xn)\0) 

is an eigenstate of the ordinary transfer matrix 

{A{fi) + D{fi))B{Xi)B{X2) . . . B{Xn)\0 ) = t(^, {Xj}) B (X^) B {X2) . . . B{Xn)\0), 

(1.11) 

and, hence, of the Hamiltonians if the parameters {Xj} satisfy the following Bethe equations 

fc=i 

These Bethe states have exactly the same property as the Bethe states for the spin ^ case 
in particular one can prove the Gaudin formula [^, 17] for their norms 

j=l k=l j^k ^ 

^'({A})a6=-^^log</.a({A}). 
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N 

and Slavnov formula |^ for scalar products of a Bethe state Y\ B{Xi:)\0) and a state 

k=l 

N 

(0 I n ^il^j) with arbitrary set of parameters: 

(oinc(M,)ni^(A.oio)^M^, ,1.14) 

j — 1 k — 1 

d 1 

T{{X, fl})ab =7TT-T"(/^6, {A}), V{{X, fl})ab = T , 

where r(//fe, {A}) is the eigenvalue of the ordinary transfer matrix ( |1.11| ). 

The thermodynamic limit for the higher spin Heisenberg model is slightly more complicated 
than in spin ^ case. The ground state of the spin s XXX model in the thermodynamic limit can 
be constructed using the "string" solutions of the Bethe equations. It was shown in Q that for 
the spin s chain the ground state contains only strings of length 2s. Such a string solution can 
be written as follows ^ 

Xl = fik + i{a- s - -), 

where Hk is real and called the string centre. The particularity of the ground state is the fact 
that it contains only strings of one particular length. 

The density of string centres in the ground state can be obtained from the following integral 
equation similar to the Lieb equation in spin ^ case: 



2s-l J °? 2s ^ 

/52.(A) + 2^ / d^iKaiX- ^)p2s{fi) + / (i//i^2s(A-//)p2.(/^) = X]^«(^ + ^(«+ 9 -^))' 

a=l_L fc=l (1.15) 



where the kernels Kj{X) are defined as 



(A + ia){X — la) 
The solution of this equation can be easily obtained 

2 cosh(7rA) 

It is a remarkable property of the spin chains that for any spin the density is the same as in the 
spin i case, but one should note here that the ground states are however very different being 
constructed from the different types of strings. 

These results for thermodynamics of the spin chains of arbitrary spin obtained by L. Takhta- 
jan and H. Babujan in |10| will be used in next sections for the calculation of the correlation 
functions. 



2 Inverse problem 

In this section we remind the solution of the inverse scattering problem for the spin s XXX 
chain. In other words we reconstruct the local spin operators in terms of the fused monodromy 
matrices. We follow the approach proposed by J.-M. Maillet and V. Terras in [^. 
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To illustrate the results of we start with the spin 1 chain. In this case the fused monodromy 
matrix has the following form: 



j A{X + i)A{X) -^{A{X + i)B{X)+B{\ + i)A{X)) B{X + i)B{X) 

r(2)(A) = 



■^{A{X + i)C{X)+ \{A{X + i)D{X) + D{X + i)A{X)+ ^(L>(A + i)B(A)+ 
+C7(A + i)^(A)) +B{X + i)C{X) + C{X + i)B{X)) +B{X + i)D{X)) 



\ C{X + i)CiX) -^iD{X + i)C{X) + C{X + i)D{X)) A{X + i)A{X) ) 

This matrix can be used to construct the operators and s^, but here we are mostly 
interested in the reconstruction of the local elementary matrices 

\ / ab 

The local operators Ej^' ^ can be expressed in terms of the corresponding monodromy matrix 
elements: 

m-l M 

E<^''^-=llr('\^,-^)Til^,JU-^) J] ^^"HCj - ^l (2-1) 

j=l j=m+l 

where t^'^^ (A) = tro T^"^^ (A) is the fused transfer matrix. 

This result can be easily generalised for an arbitrary spin s. Now the fused monodromy 
matrix is a (2s + 1) x (2s + 1) matrix and its elements are again sums of different products of 
2s fundamental monodromy matrix elements: 

^S'(^) = iL E ^^^sk,A^ + 2si-i)... T,,,,{X + i)r,,fc,(A), 

(<-^2s ^2s )2 Ji+...j2s-2»=c-l (2.2) 
fcl + ...fc2s-2s = a'-l 

where are binomial coefficients. For example the corner matrix element T^f^ [X) is just a 
product of 2s operators A{X — ki), /c = 0, . . . , 2s — 1. 

To reconstruct the elementary local operators one should again dress these monodromy 
matrix elements by the corresponding fused transfer matrices: 

m—l M 

E<.,<^^=Ylr^'^\^,-zs)4'^l,J^^-is) n (2-3) 

j=l j=m+l 

The shifts of the inhomogeneity parameters are chosen in such a way that the eigenvalue of 
the fused transfer matrix taken in the points + is on a Bethe state is: 

N . _t _■ 

r(2^)(e, - ^s)i?(Al) . . . B{X^m = J] ' B{X^) . . . B{Xm)\0). (2.4) 

This is the solution of the quantum inverse problem for the arbitrary spin Heisenberg chain. 
Now we can use these representation to calculate the correlation functions and to do it we should 
first of all understand how these complicated combination of the fundamental monodromy matrix 
elements act on Bethe states. 
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3 Finite lattice correlation functions 



As our ultimate goal is to calculate the mean values of products of local operators with respect to 
a Bethe state and, in particular to the ground state, following the same ideas as in we should 
consider the action of the local operators on Bethe states. As the Bethe states are eigenstates 
of the fused transfer matrix we need only to consider the action of the elements of the fused 
monodromy matrix. The result of this action is the "algebraic part" of the final expression for 
the correlation functions. We begin by considering the action of a single local operator on a 
Bethe state. 

First of all we remind the action of the operators A, B and D (elements of the fundamental 
monodromy matrix) on a "bra" Bethe state: 

N 

N N+l n (-^fe ~ ^a' - i) N+l 

(0 I n C{\u)A{\n+i) = «(V)^^7^^ (0| n C(Afc); (3.1) 



k=i a'=l W{Xk-K') 



k = l 
fc = l 



N 

N N+l n i^a — Afc — i) N+l 

; 1 n c{Xk) d{Xn+i) = d{Xa)'-^^ ( 1 n c{x,y, (3.2) 

n(Aa-A,) i-i 



N 



N 

N+l n (Aa - Xk — i) 

k=l 



{0\llC{Xk)B{XN+i)=YdiK 

^==1 ''=1 n (Aa - Xk) 



k^a 



N+l 

n (Aj - Aa' - i) 

N+l ^ j=i N+l 

It can be seen from these formulae that there are two different type of sums produced by the 
action of the monodromy matrix elements "j4-type" and "Z?-type", (action of the operator B 
produce for example one "j4-type" sum and one "D-type" sum). In our case the eigenvalues of 
the operators ^(A) and -D(A) in the ferromagnetic state are 

M 



One should note that in order to calculate the correlation functions we should act by the 
"strings" of operators. Consider first the action of one local operator Em°'. As it was shown 
in the previous section it can be written as a sum of ordered products of the fundamental 
monodromy matrix elements taken in the points ^ — is, ^ — is + i, . . . , ^ + is — i . The monodromy 
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matrix elements acting on the ground state produce the sums on the "D-type" indices aj and 
"A-type" indices a'j and the number of such sums is the same for all the products corresponding 
to a fixed local operator, namely a — f "D-type" sums and 2s — a' + f "A- type" sums. Introducing 
the new notations Xn+j = ^ — i{s — j + 1) we can just repeat the calculations for the spin ^ Q 
taking however into account that now d{\N+j) 7^ for j > 1: 

k=i ^ ' j=i ^ k=i _ 

TV TV ~ 

(oi n c{Xk) n B{Xk)\o) 

k=l k=l 

1 rr , , det^{{a,a'}) 

n i'^kl aj,ar=l 
k=l 

where the function H is defined as 



\TV+fcJ 



k>l 



j=l \k=l k=j+l 



a-1 TV ^ _ X _ ■ I 2s 

X n di^a,) n _^ J- n(^-. - ^^+'^ - ^) n - i ' 

j=i k=i + * \k=i k=j+i I (3.4) 

For the indices Oj, a'^ bj etc we conserve the notation of the spin ^ case: 

{61, . . . , bm} = {a2s-a+i^ . . . , a'l, ai, . . . , aa_i}. 

The determinant in the denominator is the Gaudin determinant and the matrix in the numerator 
is also the Gaudin matrix with some replaced columns. We will consider this "analytic part" in 
general later on in this section, and for the ground state for the spin 1 case in the last sections. 
One can easily obtain representations for the replaced columns from the scalar product formula. 
Here we will consider the "algebraic part" of the expression (3.4). Taking into account values of 
the parameters Xjsi+k we obtain: 

j=i \k=i k=j+i ) 

One can see here that the product corresponding to the "D-type" parameters is the same for 
any element of the sum in ( |2.2| ) . Moreover from this result one can see that aj > N gives a non 
zero contribution only if > + 1 and there is Uk = aj — 1, k < j. It leads to the conclusion 
that such a term should contain d(^ — is) which is zero. It means that as in spin ^ case the 
summations over aj should be taken only from 1 to A^. 

The product corresponding to the "A-type" parameters is not the same for all the terms but 
taking the sum in (2^) and symmetrizing over the permutations of the "I?-type" and "^-type" 
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parameters separately one can simplify it and obtain finally: 

N N 

(0| UC{\k)ET UB{\k)\0) 

//-, \ k=l k=l 

= — 

{0\UC{Xk) UB{Xk)\0) 

k=l k=l 

-( \- \- Txa a . det*({aX}) 

n i^k\ «j=i 
fe=l 



a— 1 , , 2s— a+1 \ \ a— 12s— a+1 



Hf.,..;)(w>=n (/!!;pVi n TXT^xTPTin n , 

fc,i=i ^'^"fc '^"i'' ^-^ fe,i=i '^'^"S ^O'V ^ k=i 1=1 

fc>i fe>! 
a-1 /2s-l \ 2s-a+l /2s-l \ 

^WiWiK-^+^i^-m n \ X{^^a'.-i+i{s-k+i))\ 

3=1 \k=l J 3=1 \k=l J (3.7) 

One can see from this representation that terms with a'^ > N produce non zero contribiitions if 
a'j = N + 2s or if there is = a'j + 1. It means that the operators — i{s — k)) which appear 
in the scalar product after the action of the local operators should form a "substring" without 
holes starting from ^ — is. For example states like 

{0\C{^-is)C{^-i{s-l))...C{^-i{s-k)) n 

b<N 
^ J' J 

produce non zero contributions to the correlation functions but the contribution of such states 
as 

{0\C{^-ts)C{^-i{s-2))l[CiXh) or ( - i(s - 1)) [] C(Ab) 

b<N b<N 

b^a^.a'. b^aA.a'. 
3 3 1 ] 

are zero. This property is rather important as the matrix appearing in the scalar product is 
much simpler in this case. 

Consider the determinant det ^ appearing in (3.4) from the scalar product 

N N 

( |C(e - is)C{i - i{s -I))... C{i -i{s-k + 1)) n C{\) n ^(Aa)| 0) 

b=k+l a=l 

After extracting of the normalisation coefficients we obtain the following matrix: 
^ab =<6, h > k, 

*"'"(A„-^ + is)(A„-^ + z(5-l))' 

^{Xa-^ + i{s-b+l)){Xa-^ + i{s-b)) + 

^^\Xa-^ + i{s-b+2)){Xb-^ + i{s-b + i)y ^<^^^- 
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where 

The columns with 1 < 6 < A; can be considered as sum's of two columns and the second one 
does not contribute to the determinant as it is always linearly dependent on the columns with 
a' < a. Thus the matrix ^ can be replaced by ^ 

*a& =<6, h > k, 

*a& =7T T— 7 , . A. . ——, TTT = PbiXa - 0, b< k. 

{Xa-^ + i{s-b + l))(Aa -^ + i{s- b)) 

A similar calculations can be done also for m-point functions leading to the following repre- 
sentation: 

N m , N 

(oinc^(A.) n^r ' n^(A.)io) 

k=l 1=1 fc=l 



(oi nc(Afe) n^(Afc)io; 



fe=l k=l 



=(_l)EK-i) 



n c^r'cT 



/2s-l m 2s 2r-l / n 

n ^'^fe! n n n 

\fc=l / j,fc=i r=l n=l ^ ^ 

3>fc 

X i: E Hf,.,,,({A})^5i|i|^. (3.8) 

Here we introduced some new notations. In every site I there is a local operator which 
produces sums over ai — 1 "D-type" indices aji and 2s + 1 — "A-type" indices a^^, and we 
define Ajv+2si+fe = 6 — '^(s — A;+l). We obtain the algebraic part: 

m 2s-l fa-i-l 2s+l-a[ \ 

H?;s.j({A})=nn n(^«.<-^'+^(^-^)) n (A„;,-6+^(^-^+l)) u 

/=! fc=l V-?=l / 

X n ^on})Gtn{{^'jU «.n})Gf/({a;,, a^-J) 

n<i (3.9) 
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Where the "two sites" contributions for two different sites I and n are 

a; — 1 an — I 

n (^s, -^n- is) n (^afc„ - 6 + ^s) 

G^"" =— — , (3.10) 

ai-la„-l ^ ' 

n n {Kii - Kn - i) 
j=i k=i 

— 1 2s—an+l 
U {Xa,, - - is) U {Xa',^-Cl-i{s-l)) 

^^^^ =~ a,-12.-<+l " ' (^-11) 

n n {Xa,-Xa',„-^) 
j=l k=l 

2s-q;+1 a„-l 

n (Aa^, - + i(s + 1)) n {Xa,^-^l+is) 

=^=^^ wTl-^ ^^^^ ' (3-12) 

n n (Aa.„ - A„^.^ - 
j=i fc=i 

n (A„.., n (A„.^-Ci-i(5-i)) 

Gi^ =^ -, '-^ , (3.13) 

n n (A4„-Aa— 

j=l fe=l 

and the diagonal terms are given by 

^« ^« - 11 (A„.^-A„,J2 + i 11 (at^aTT + T' ^ ^ ^ 

Q,,_l2s-aJ+l 

We obtained here the algebraic part of the expression for a m-point correlation function for 
an arbitrary spin Heisenberg chain. Considering the determinant part we can easily see using 
the same arguments as in one-point case that the main proposition about "the substrings" (the 
terms in the svuti (3.8) with a'ji > N produce non zero contributions if a'ji = N + 2sk < / or if 
there is a'^^ = a'ji + 1 with k < j) is also valid in the m point case. It means that the determinants 
are always simple and contain only Gaudin columns and columns Pa(A(, — ^k)- More precisely: 
every term of the sum (3.8) contain a determinant which is obtained from the scalar product 

M ki N N 

1=1 j=l 6=A;i+---+fe,„+l a=l 

where ki = a'- , — 1, the "substring end" a'- , being the minimal a'-, > N in the corresponding 
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term. The corresponding matrix appearing in the sum (3.8) has the following columns 



1=1 

^ l~l I 



{Xa-^l+i{s-j + l)){Xa-Cl+i{s-j)) 
/-I 



where j = b — ^ k^. 



r=l 



Thus we have a representation for the correlation functions for a finite arbitrary spin Heisen- 
berg chain. Being in some sense very similar to their spin | counterparts these representations 
are rather complicated for big spins. For this reason in the next sections we will consider only 
the first generalisation of the results of Q which is the spin 1 chain. 

The next step of our approach is the thermodynamic limit for the ground state. On this 
stage the main difference with the spin ^ case appears as the ground state is constructed of 
bound states (2s-strings). It produces some new difficulties which will be considered in the next 
section using the simplest example of the one point functions. 



4 One point functions 

To illustrate the last stage of the calculation, i.e. the introduction of the string solution of the 
Bethe equations for the ground state we begin with the simplest example, namely with the one 
point functions (corresponding to the diagonal elementary matrices E^^) 

Consider the simplest correlation function (1-point emptiness formation probability) of the 
spin 1 XXX chain in the homogeneous case: 

/3(1) = {vac\{Tf^)-\-i)D{Q)D{-i)\vac) (4.1) 

Using the action of the operators D on the vacuum and the scalar product formula we easily 
obtain (there is no difference with the case spin 1/2) for the finite chain: 

a b^a 

Where is the Gaudin matrix : 

= {MK{Xj) - J2 K{\j - \i))5jk + K{\, - \k) (4.3) 

with 

2 

^(^) = (A + z)(A-z) 
and the matrix ^ is obtained from the scalar product: 

^bk = . , =P'-{W- 
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Hence one can again divide one matrix by another and reduce it to a 2 x 2 matrix: 

where (p^ are solution of the following systems of linear equations: 

(MKiXk) - - Xk))<Pt + Yl ^(^J- - ^i^^^^f = P'±(^k) (4.4) 

j j 

Thus the 1-point function (|4.1| ) is given by 

/3(i) = iyy ^ det f 5 (4-5) 

In the thermodynamic hmit the ground state of the spin 1 XXX chain is built of 2-strings 

M->CO: A2A:-1 ^ ^fc + -, A2A:^AiA:--, Im(/ifc) = 

To obtain the equations for the analytic part in the thermodynamic limit one should take into 
account the finite size corrections to this string picture as some terms in ( |4.4| ) become singular. 
To analyse the excited states with finite energy of the XXX spin ^ one usually considers first 
the string limit and only then the thermodynamic limit, as these corrections are exponentially 
small. However for the ground state of the spin 1 XXX model one cannot use this method in 
a rigorous wayQ as the corrections to the string picture calculated in |19, 20, 21, ^ are of the 
order 

Ot]^ 1 

A2fc-i - A2fc - i = 2i— + o(— ), a>0 

The correction a is always positive and it makes possible to rewrite the system of linear equations 
( [4.4| ) as integral equations in the thermodynamic limit with a special choice of the integration 
contours near the singular point: 

oo oo 

j dXKip- X)^fiX)+ j dXK{n- X + i + iO)ipf{X) =p'^ip + ^), 



-oo — oo 

oo oo 



iffifi) + / dXKif, - X)^fiX) + I dXK{i^-X-i- iO)ipfiX) =p'±(/x - '-), (4.6) 

where ff^fik) = Mp{fj.k)(j)^j^_i and Lp^il^k) = Mp{fj.k)'P^k- used here the integral equation 
for density of strings in the ground state 

P(A) ^ 



2 cosh(7rA) 



^this method leads to the same result for the correlation functions as one described below but in a more 
complicated way 
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The solution of this system can be easily obtained 



(/.+ (A)=p(A), v^^(A) = 0, 

ipi{\) = 0, ^^{\)=p{X). (4.7) 

Now we can substitute this result to the expression for the one-point function ([4.5|), replacing 
sums by integral avoiding the singular point in the same way as in the integral equation (taking 
into account the sign of the finite size corrections) (4.6). 



oo 



i f f ... 1 / (A + !)(/. -I) (A-i)(^ + i) 

cosh(7rA) cosh(7r/i) I A — /i + iO A — /x — 2i 



/3(1) I dX I df. 



(4.8) 



One of these two integrals can be calculated as only the pole in the point fi = X contributes. 
We obtain finally: 

— oo 

Of course this result can be obtained directly from the symmetry of the model, but this cal- 
culation illustrate well how to deal with strings in our method and it can be useful not only 
for the more general case of m point functions which will be considered in the next section but 
also for the computation of more general correlation functions, depending, for example, on the 
temperature. 

Two other one point functions can be calculated in a rather similar way but here we should 
consider also the "A-type" sums which contain more terms than "D-type" sums considered in 
the previous example. We will show, using the simplest example of the one point functions 
/2(1) and /i(l) that this problem can be solved exactly as in the spin ^ case by moving the 
corresponding contour of integration. 

Consider the function 

/2(1) = ^{vac \T2\-i) (^A{0)D{-i) + D{0)A{-i) + C{0)B{-i) + 5(0)C(-i)) | vac) 
We easily obtain a finite lattice representation for this function (3.6) 

One should mention that the sum over index a is a L>-type sum and the sum over index 6 is a 
^-type sum and it contains one additional term (only one because of the substring limitation). 
As in the previous case we can rewrite these sums as integrals in the thermodynamic limit: 



/2(1) = -'- I dX I df. 



1 /(A + i)(^ + i) (A-i)(/z + 



cosh(7rA) cosh(7r^) \ X — + iO A — /x 



-00 
00 




+ / dX ^. ^, ^ ^, (4.11) 
' cosh(7rA) Aa - 3§ 
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One should note that function p(A) = 2 cosh (tt A) same as in the spin i case and has a pole 



at A = — I and its residue is 



27ri Res/9(A)|^^ 



It means in particular that shifting the contour of integration on the variable jjl to the line 

parallel to the real axis with Im(yu) = —1 one crosses the pole of p(A) for the second double 
integral in (4.11) (for the first double integrals one can move the contour without crossing any 
poles) and the contribution of the pole is exactly the single integral in (4.11). Finally we get 

00 00 / . . . . \ 

^ f- 1 AA + |)(Ai-|) (A-|)(M+|)\ 



/2(1) J d\ J d^i 



cosh(7rA) cosh(7r/i) I X — /i + i X — fi — i 

—00 —00 ^ 

Here once again one can reduce this expression to a single integral: 



00 



TT 



i - 2A2 1 



— 00 

Considering the last one point function /i(l) one can obtain in a similar way that one should 
move both contours to obtain the same double integral representation as for /3(1): 

—00 

In this section we considered the simplest examples of the correlation functions, but however 
this simple example illustrates quite well the basic properties of the thermodynamic limits for 
the ground state constructed of the 2-strings. We have also shown that the ^-type sums should 
be replaced by the integrations over shifted contours (as in spin ^ case). 

In the last section we show how this technique can be used for the general m-point functions 
for the spin 1 Heisenberg XXX chain. The proofs of them in general are absolutely equivalent 
to the calculations in this section, but contain some very cumbersome formulae which we omit 
in general for the intermediate steps. Also this method can be used for arbitrary spin and we 
give in the end one of the possible multiple integral representations for the correlation functions 
of the higher spin chains. 



5 Correlation functions 

In this section we generalise the results of the previous section for the general m-point equal-time 
correlation functions of the XXX chain spin 1 in the thermodynamic limit. 

We calculate the following correlation functions or, more precisely, the elementary blocks 
which permit to construct any m point correlation function: 



f{a,a'}{m) 



in „ . 



(5.1) 
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where E^^' ^ are elementary local 3x3 matrices -E;'^'" = Si^a'^k,a and \ipg) is the ground state 
of the model (in the spin 1 case the number of quasiparticle N in the ground state is equal to 
the number of sites M). 

For a finite spin 1 chain we obtained (3.8) that this correlation function can be represented 
as multiple sums: 

/ 1 1 \ ^ 

= (-l)S('^J-^) '-^ ^ H{„,„,}({A})det2„.S({a,a'}). 

n (Ci - ^k^m - 6)' + 1) {aj,o^} (5.2) 

i,fc=i 
j>k 

In the thermodynamic limit the solution of the Bethe equations corresponding to the ground 
state consists of 2-strings distributed with the following density: 

1 ^ 

Ptot(A) = J^^p{>^-Cn), 

n=l 

As in the spin ^ case it is convenient to introduce a set of indices bk which is the set of aji,a'ji 
ordered in a special way: 

{h,---, b2m} = {{a'jm}^ • • • > Wjl}: {aji}, {ajm}}, 

with local subsets {aji}, I < j < ai — 1 and {OjJ, 1 < J < 3 — a'l, (these subsets can contain 
one, two or no elements). 

Let us now consider the determinant part of the expression (3.8). In the thermodynamic 
limit we can divide the matrix in the numerator by the Gaudin matrix, or more precisely we can 
calculate det($'~^^'({a,a'})). This determinant can be written as a determinant of a 2m x 2m 
matrix and the matrix elements of this matrix are given by the inhomogeneous version the 
integral equations (4.6) with p±{ii — ± |) in the r.h.s. Due to the translation invariance of 
this equations we obtain just the same solution with a shift 2!/^ ~ Finally we obtain: 

detM($'"'*({a,a'})) 
Sjk{{a,a'}) 

Sj2k-ii{a,a}) 
Sj2k{{a,a'}) 

Now all the sums over bi from 1 to M can be replaced by integrals taking into account 
account that we obtain a sum of two integrals obtained from bi = 2j, j = 1, . . . M/2 and 
bi = 2j — 1, j = 1, . . . M/2 dealing with the contours near the singularities of the algebraic 
part in the same way as in the previous section. Replacing sums by integrals we use the following 



:det2m'5({a,a'}), 

■■ - Sbj-M,k, bj > M, 

1 + p(A6. - Cfc - i) 



2 Ptot{\ - I) 
1 _(_!)&.■ p(A,. - 4 + i 
2 p,,,{X, + i) 



-, bj < M, 



■, bj < M. 



(5.3) 
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rules and notations: 

M/2 

E - 

i,;=2j-l 



M/2 

E 

<'fcr-=2j-l 
J = l 

M/2 

E 



-oo + l+iO 
oo+l+iO 



C^i^fcr Ptot{l^'kr 



M/2 

E 

bi=2j 

i=i 
M/2 

E 

3 = 1 

M/2 

E 

J = l 



C^A; Ptot(Az + 



C^Z^fcr Ptoti^kr + 2) 



-oo-|-iO 



oo-|-iO 



-oo-|-iO 



In this notation the set {A;} is the same set as {ffer, ^kr) but ordered in a special way: 

{Al, . . . , \2m} = {Wkm}^ ■ ■ ■ > Wkl}^ {z^fcl}, ■ ■ ■ , {ffcm}}- 

We will use also the following parameters: 



ei =2, 



Q = -2' 



if A; = 1/kr (D - type), 
if A; = ffc^ {A - type), 



we can also associate to every A; the corresponding site number r; if A; = v^ri 01c Xi = u'^^^ . This 
notation are very useful to simplify our formulae. 

Replacing the sums by two integrals one should mention that for different contours for \j we 
obtain different functions in the determinant, namely —^k ~ §) for the upper contour and 

y?^(Aj — ^jfc + I) for the lower one. To simplify the formulae we introduce the following function: 



<A±(A)=(^f(A--), Im(A)> 



</.±(A)=^±(A + ^), Im(A)<-^ 



1 

2' 
1 



(5.4) 



(Note that we can deal with this function as with an analytic function if the integration contours 
do not cross the line Im(A) = —\)- 

We should now analyse the terms with a^,^ > M. As in the spin ^ case they can be written as 
integrals around the corresponding poles of the determinant part in the points Xj = ^fc (for the 
upper contour), Xj = ^jt — i (for the lower contour). More precisely we can replace the complete 
sum over a'j^^ > M by the following integral: 



E 



E 




dT'kr PU>t{l^'kr + 2^ 
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where F; and T'^ are small contours around the points ^; and — i. It means that as in the 
spin ^ case the eontribution of the poles can be absorbed into the integrals by moving all the 
contours for z^^^ (A-type variables) down by i (to avoid crossing of some additional poles one 
should move first all the lower contours and than the upper ones). 

Now we can write a multiple integral representation for the correlation functions of the spin 
1 XXX chain: 



f{a,a'}im) = i"' 



1 



'=1 - X 

m ^ 
3>k 



X 



j + j |dAi..| J + J |ciA2mH{„,„,}({AJ)det2m5({A}). ^ 



The algebraic part for spin 1 can be written as: 



2m m 



H{a,a'}({Ai})=nn(^'-^*^ 



2m I ri — 1 m 

n n (A; - 6 - 2i£/) n {k-^k + ^iei) 

1=1 \ k=l k=ri+l J 



l>n \ J 

and the 2m x 2m matrix S is defined as 

>S,-,2fe=</>+(A,--6)- (5.7) 

This is one of many possible forms of the results which can be easily generalised to the higher 
spin cases. One should note that the correlation functions arc written once again as a multiple 
integral (or a sum of multiple integrals) and the integrals are taken over the solution of Bethe 
equations for the ground state (2-strings here). The expression under the integral can be once 
again separated into two distinctive parts: one defined by the choice of local operators (algebraic 
part) and one defined by the ground state (determinant). 

This result can be also rewritten in many different forms. First of all as for the one-point 
functions one can reduce the number of integration and obtain the result only as m integrals 
(this is a particularity of the spin 1 case, it cannot be done for the higher spins). First of all one 
should note that Sj^2k-i = for upper integration contours and Sj^2k = for lower ones (4.7). 
It means that the determinant has a block diagonal structure and the integrals can be rewritten 
as a sum of integrals with m integrals over the upper contours and m integrals over lower ones: 

f{a,a'}{m) = — ^ ^- J2 (-1)'"^ / ^/^l--- / '^^-^ 

n (0 - ^kM^j - + 1) {A}=Mu{.'} _j _j 

j,k=l ^ ^ 

j>k 

)+i+iO oo+i+iO 

I dfi[--- I d/x;„H{,,,,}({AJ)det^W({/x+|})W(V-|}). (5.8) 
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where the sum is taken over all possible partitions of the set {A} with 2m elements into two 
subsets with ra elements and [a] is the sign of the following permutation 



Cr({Al, . . . , \2m}) = {/^l, • • • ,/im,/^l, • • • 

and det>V({;u'}) is the spin ^ determinant: 

y^jk = p{fJ-j - Ck)- 

this representation is also convenient as it contains only m x m matrices and only meromorphic 
functions. Now we can move the contours in this integral to obtain that only poles in Xj — 
Afc + i{ej + ek) contribute, which gives a representation of the result as a sum over all possible 
splitting of the set of 2m variables A/ into m pairs which form strings A/ = z^j — |, A^ = vj + | 
but every term contains only m integrals over the string centres vj: 



f{a,a'}{m) =(27r)" 



n cr'c^' 



1=1 



— j dur-- j du^ det^>V({z^}) x 



j>k 



X Yl G|„,„,|({z/},{{Z,Z;.},j = l...m}). (5.9) 

{l,2,...,2m}=UJL^{ljir} 

The algebraic part here is however much more complicated than in the spin ^ case. It is 
written as a sum of (2m— 1)!! = 1 • 3 •. . . • (2m — 1) terms, where every term is the corresponding 
residue of the general algebraic part 

G{a,a'}{M,{{lj^'j}^j = l...m}) = 

= Res|,^^^,^, +,(,^^+,^, . . . Res|,^^^,^^+,(,^^^,^,^). H|,,,,}({AJ). (5.10) 

This is another possible form for the final result. It is rather particular as the correlation 
functions are represented as m-integrals as in the spin ^ case and the analytic part is just the 
square of the analytic part for spin i. It is also important to note that the homogeneous limit 
of this expression can be obtained exactly as in the spin ^ case. 

d<wi{,}) ^ ^^^^ w^°-m), wfr = (,'.!!T-i PiH) 

3>k 

Once again the analytic part will be just the square of the corresponding spin ^ algebraic part. 

Thus we obtained several equivalent expressions for the correlation functions of the spin 
1 Heisenberg XXX chain. It is interesting to note that this result looks quite different in 



comparison to the corresponding results obtained in [13, |lj] (it is quite clear that for the same 
quantity one can write many different integral representations and sometimes it is rather difficult 
to prove that they are equivalent). The first (5.8) result has a form which is much simpler, even if 
now we have 2m integrals instead of m. Moreover similar results can be obtained for higher spins 
(3.8) in the thermodynamic limit, with more complicated contours (in some sense we integrate 
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always over the strings), determinant of a 2ms x 2ms matrix and corresponding general algebraic 
part: 



f{a,a'}im,s) 



1 

™ 1 ^' 1 N 2 

n c^r'cT 
1=1 



/2s- 1 \ m 2s 2r-l . 

n ^'ki) n n n fe-a-^ir-n 

\fc=l / j.fe=i r=l 71=1 ^ 

j>k 



j>k 

X 

2s-strings 2s-strings 

Here the intergals are taken over the strings with corresponding little shifts: 

oo-(i+i0)(2s-fc-|) 



(5.12) 



I dX /(A) = f2 I dX /(A), 

2s-; ''^^ 

the algebraic part is 



2s-strmgs "'"-oo-(i+i0)(2s-fc-i) 



2s— -1 2sm m 



Hit'}({A/}) = I n n n(^' - - - 



p=i 1=1 k=i 

2sm I ri — 1 m ' 

n n (A/ - - 2isei) n (A^ -Ck + ^iseO 

1=1 \ k=l k=ri+l 

X ^ ^ (5.13) 

n ( Az - An - i(ez + en)^ 

l>n \ 



and the 2sm x 2sm matrix 5^*^ is defined as 



Sj^2s{k-i)+i =p{^j -ik + 'isi-il + -), - 2s + / - 1 < Im(Aj) < -2s + I 

=0, otherwise. (5-14) 

This representation can be proved in a very similar way to the spin 1 case (5.8). 

We discussed in this paper the correlation functions of the higher spin XXX chains and we 
have shown that the even if the ground state contains bound states the correlation functions 
can be calculated. We hope to use these result to calculate the mean values of local operators 
with respect to any excited state for the spin ^ Heisenberg chains and, hence, to obtain a 
representation for the finite temperature correlation functions. 
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